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The  angular  momentum  propagated  by  a  beam  of  radiation  has  two  contributions:  spin  angular 
momentum  (SAM)  and  orbital  angular  momentum  (OAM).  SAM  corresponds  to  wave  polarisation, 
while  OAM-carrying  beams  are  characterized  by  a  phase  which  is  a  function  of  azimuth.  We  demon¬ 
strate  experimentally  that  radio  beams  propagating  OAM  can  be  generated  and  coherently  detected 
using  ordinary  electric  dipole  antennas.  The  results  presented  here  could  pave  the  way  for  novel 
radio  OAM  applications  in  technology  and  science,  including  radio  communication,  passive  remote 
sensing,  and  new  types  of  active  (continuous  or  pulsed  transmission)  electromagnetic  measurements. 

PACS  numbers:  41.20.-q,  41.20.Jb,  07.57.Hm,  07.57.Kp, 


As  early  as  1909  the  wave  motions  associated  with  a  re¬ 
volving  shaft  were  studied  by  Poynting,  who  suggested, 
by  analogy,  that  circularly-polarised  light  should  carry 
angular  momentum  (OAM)  [?  ].  In  1936  Beth  reported 
on  an  optical  experiment  [1]  in  which  he  verified  Poynt- 
ing’s  prediction.  Since  then  others  have  studied  electro¬ 
magnetic  OAM  [2?  ,  3],  however  it  is  only  recently  that 
electromagnetic  OAM  has  been  extensively  studied  and 
utilized,  both  theoretically  [4?  -6]  and  experimentally,  at 
visible  [7?  ,  8],  millimeter  and  microwave  [9?  12],  and 

radio  wavelengths  [13?  ].  OAM  has  also  been  studied  in 
electron  beams  [14]. 

In  addition  to  energy  and  linear  momentum,  electro¬ 
magnetic  waves  can  propagate  angular  momentum  to  in¬ 
finity.  The  total  electromagnetic  angular  momentum  has 
mode  numbers  j  =  s+l ,  where  s  and  l  are  mode  numbers, 
respectively,  of  the  electromagnetic  spin  and  orbital  an¬ 
gular  momentum[15].  Classically,  spin  angular  momen¬ 
tum  manifests  itself  as  wave  polarisation,  where  mode 
numbers  s  =  ±1  correspond  to  right-  and  left-hand  cir¬ 
cularly  polarised  modes  and  s  =  0  to  linearly-polarised 
modes.  Similarly,  OAM-carrying  beams  are  character¬ 
ized  by  a  phase  which  is  a  function  of  azimuth,  and  OAM 
mode  numbers  l  are  classically  manifested  by  a  change  in 
phase  of  l  x  360°  around  any  arbitrary  circle  centred  on 
the  beam  axis  [7,  11]. 

We  extend  the  radio  angular  momentum  technique 
[2,  5,  6,  9,  10,  12,  13,  16?  T9]  by  actively  generat¬ 

ing  OAM  radio  beams  having  a  variety  of  mode  num¬ 
bers  using  a  simple  antenna  array  [5],  and  by  using  a 
phase-coherent  technique  to  detect  and  verify  the  trans¬ 
mitted  modes  [19].  The  results  presented  here  demon¬ 
strate  that  radio  beams  propagating  OAM  can  be  gen¬ 


erated  [5,  6]  and  coherently  detected  [19]  using  ordinary 
electric  dipole  antennas.  This  could  potentially  pave  the 
way  for  novel  radio  OAM  applications  in  technology  and 
science,  including  radio  communication,  radio  and  radar 
remote  sensing,  ionospheric  radio  diagnostics,  and  radio 
astronomy.  Note  that,  using  current  technology,  it  is  only 
at  radio  frequencies  that  coherent  measurements  of  elec¬ 
tromagnetic  fields,  i.e.  measurement  of  both  amplitude 
and  phase,  may  be  performed. 

An  antenna  array  of  N  identical  sources  at  angular 
frequency  to  and  equal  amplitudes  has  the  array  factor 

N 

V  ^  exp[— i(fc  •  xn  -  <j)n)\,  (1) 

n—l 


where  k  is  the  wave  vector,  xn  is  the  position,  <f>n  is  the 
phase  of  the  nth  emitter,  and  i  represents  \AM[6,  20- 
22].  When  the  emitters  are  electric  dipoles  with  dipole 
moments  d,  the  electromagnetic  energy  density  u  = 
£o(\E\2  +  c2\B\2)/2,  where  E  and  B  are  the  electric  and 
magnetic  fields,  becomes 


k2\k  x  d  2|T|2 
e0(47rr)2 


+  0(r-3) , 


(2) 


where  £o  is  the  vacuum  permittivity  and  r  =  |a?|  is  the 
radial  distance  from  the  centre  of  the  array,  where  x  is  the 
position  vector.  The  total  angular  momentum  density  is 
defined  as  h  =  x  x  g,  where  g  =  £oR,e{E  x  B*}  is  the 
linear  momentum  density.  In  our  case  we  have  linearly- 
polarised  dipoles,  i.e.  d  x  d*  =  0  and  s  =  0,  so  that 
the  total  angular  momentum  j  =  l  +  s  =  l.  Hence,  only 
the  OAM  contributes  to  the  total  angular  momentum 
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density, 


h  = 


u  m>  *M> 

u  |T|2 


+  0(r"3), 


(3) 


where  L  =  —i(x  x  V)  is  the  OAM  operator.  Since  the 
OAM  operator  in  Eq.  3  only  operates  on  the  array  factor 
'T,  the  only  contribution  to  the  angular  momentum  of 
the  fields  will  be  from  the  phasing  of  the  elements  and 
the  geometry  of  the  array,  not  the  individual  antenna 
elements,  so  long  as  the  elements  are  much  smaller  than 
the  size  of  the  array. 

We  have  used  a  circular  array  of  N  =  6  elements  placed 
in  the  xy  plane  with  emitters  distributed  equidistantly 
around  the  perimeter  of  a  circle  of  radius  a  as  shown  in 
fig.  1,  and  phased  such  that  </>n  =  2Tiln/N,  with  l  an 
integer.  The  array  factor  becomes 

=  N(—i)1  exp(ilip)Ji(kasmd),  (4) 

where,  6  and  ip  are  the  spherical  polar  and  azimuthal 
angles,  respectively,  and  J;  is  the  Bessel  function.  The 
array  factor  contains  the  pliasor  exp(i^)  which  is  a 
characteristic  of  OAM  beams[7].  The  phase  factor  varies 
around  the  beam  axis  and  the  OAM  number,  l,  corre¬ 
sponds  to  a  Fourier  component.  Thus,  by  measuring  the 
phase  of  a  single  field  component,  the  OAM  modes  can 
be  separated  by  a  spatial  Fourier  transform  about  the  2 
axis  [23]. 

An  array  of  N  =  6  folded  half- wavelength  dipole  anten¬ 
nas  was  constructed,  as  shown  in  the  inset  in  fig.  1.  The 
antennas  were  fed  with  equal  amplitudes  and  with  phase 
shifts  of  5<f>  =  0°,  ±60°,  and  ±120°  between  consecutive 
elements,  such  that  the  array  generated  beams  carrying 
OAM  modes  l  =  0,  ±1,  and  ±2,  respectively.  The  radia¬ 
tion  patterns  of  the  simulated  beams  are  shown  in  fig.  2. 
In  the  simulations,  ordinary  dipoles  were  used  instead  of 
folded  dipoles  for  simplicity. 

As  can  be  seen  in  the  phase  plots  in  the  upper  panels 
of  fig.  3,  when  looking  up  the  beam  towards  the  trans¬ 
mitting  array  (i.e.  towards  the  page  in  fig.  3),  the  lines 
of  constant  phase  spiral  out  from  the  beam  axis  coun¬ 
terclockwise  for  l  >  0  (right-handed)  and  clockwise  (left- 
handed)  for  l  <  0.  One  complete  turn  at  constant  radius 
around  the  beam  axis  changes  the  phase  of  the  field  com¬ 
ponents  by  l  x  360°.  Accordingly,  for  the  measurements 
presented  here,  the  phase  of  Ey,  which  is  the  strongest 
electric  field  component,  was  measured  at  a  constant  ra¬ 
dius  from  the  beam  axis  [19]. 

The  measurements  took  place  in  an  anechoic  chamber 
at  a  frequency  of  2.383  GHz  (corresponding  to  a  wave¬ 
length,  A,  of  0.126  m)  and  at  a  distance  of  z  =  4.82 
m  (38.3  A)  from  the  transmitting  array  and  hence  well 
into  the  far-held  region.  The  Ey  held  was  measured  with 
a  half-wave  dipole  antenna  placed  at  the  tip  of  a  one- 
meter-long  rotating  arm,  which  was  mounted  on  a  pillar 


FIG.  1.  (Colour  on-line)  The  experimental  setup.  The  six- 
element  transmitting  antenna  array  is  at  the  lower  left,  at 
the  origin  of  the  coordinate  system,  with  a  front  view  of  the 
array  shown  in  the  inset  at  top  left.  The  receiving  antenna 
is  located  to  the  right,  at  the  end  of  a  rotating  arm.  The 
arm  was  connected  to  the  turntable,  the  pivot  of  which  was 
located  on  the  z  axis,  and  which  was  supported  by  a  ver¬ 
tical  pillar  (blue)  connected  to  the  floor.  The  transmitting 
antennas  were  half-wavelength  folded  dipoles  directed  along 
the  vertical,  parallel  to  the  y  axis,  positioned  equidistantly 
around  a  circle  of  radius  a  =  0.75  A  (where  A  is  the  wave¬ 
length)  and  placed  parallel  to,  and  at  a  distance  of  0.013  m 
(0.1  A)  in  front  of,  a  conducting  plane.  The  xy  plane,  and 
the  dashed  circle  of  the  transmitting  array  (see  inset),  were 
parallel  to  a  chamber  wall,  and  the  z  axis  was  perpendicular 
to  the  same  wall.  By  feeding  the  transmitting  antennas  with 
equal  amplitudes  and  with  a  phase  difference  between  consec¬ 
utive  elements  of  0°,  ±60°,  and  ±120°,  OAM  modes  l  =  0, 
l  =  ±1,  and  l  =  ±2,  respectively,  were  generated.  The  centre 
axis  of  the  horizontally-transmitted  beam  was  aligned  with 
the  rotation  axis  of  the  rotating  arm.  The  receiving  antenna 
was  attached  to  the  arm,  so  that  it  rotated  around  the  beam 
axis  at  a  radius  of  1.00  meter,  thus  avoiding  the  null  at  the 
centre  of  the  OAM  beams.  The  measurement  plane  was  lo¬ 
cated  4.82  m  (38.3  A)  from  the  transmitting  array;  the  angle 
6  was  11.7°.  During  the  measurements,  the  arm  was  rotated 
in  a  vertical  plane,  parallel  to  the  xy  plane.  The  receiving  an¬ 
tenna  was  vertically-aligned  and  measured  the  y  component 
of  the  transmitted  electric  field. 


attached  to  the  floor,  see  fig.  1.  The  amplitude  and  phase 
of  Ey  were  measured  every  4°,  for  l  =  0,  ±1,  and  ±2. 
Repeated  measurements  for  l  =  1  showed  that  the  beam 
amplitude  and  phase  were  stable. 

For  l  =  0  the  measured  data  show  a  sinusoidally- 
oscillating,  rather  than  constant,  phase  around  the  mea¬ 
surement  circle.  This  oscillation  was  found  to  be  pro¬ 
duced  by  a  small  misalignment  between  the  centre  of 
the  transmitted  beam  and  the  centre  of  the  measure¬ 
ment  circle.  This  misalignment  also  produced  a  periodic 
phase  oscillation  in  the  l  =  ±1  and  ±2  data.  The  l  =  0 
measurement  was  used  to  correct  the  phase  data  for  all 
modes.  The  ripples  in  the  phase  plots  in  the  lower  panels 
of  fig.  3  indicate  the  effect  of  spatial  reflections  of  the 
signal  and  that  the  transmission  was  not  a  pure  OAM 
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FIG.  2.  (Colour  on-line)  Simulated  radiation  patterns  for  the  l  =  0,  ±1,  and  ±2  radio  beams.  For  l  ^  0  the  beams  exhibit 
the  characteristic  null  along  the  beam  axis  ( 9  =  0°)[5,  6,  19,  24].  Simulations,  analytical  calculations,  and  measurements  show 
that  the  1  =  0  and  ±1  beams  have  maxima  at  6  =  0°  and  6  =  22°,  respectively.  For  l  =  ±2,  the  angle  of  the  beam  varied 
significantly  around  the  beam  axis.  This  is  because  l  =  ±2  is  much  closer  to  the  theoretical  limit  in  the  inequality  |Z|  <  N/2 
for  unambiguous  transmission  of  OAM  modes;  the  OAM  beam  will  degrade  as  the  limit  is  approached.  Note  also  that  the 
amplitude  patterns  for  l  =  ±2  exhibit  ripples.  Small  ripples  were  also  present  for  l  =  ±1  but  are  not  visible  in  the  figure. 
Ripples  generally  appear  when  |Z|  approaches  the  upper  limit  of  N/2.  The  color  scale  and  contour  lines  indicate  power  density 
for  each  beam  relative  to  the  maximum  of  the  l  =  0  beam;  the  interval  between  contour  lines  is  3.2  dB. 


FIG.  3.  (Colour  on-line)  Simulated  and  measured  phase  distributions  in  the  OAM  beam.  In  the  upper  panels  the  positive  a 
axis  is  out  of  the  paper,  in  the  lower  panels  the  positive  z  axis  is  oriented  upwards.  The  upper  panels  show  the  simulated  OAM 
beam  phase  distribution  for  the  y  component  of  the  electric  field  at  the  measurement  plane.  From  left  to  right,  the  panels  show 
the  phase  for  l  =  —2,  —1,  0,  1,  and  2,  as  seen  looking  along  the  beam  axis  toward  the  transmitting  array,  i.e.  in  the  negative  a 
direction.  The  white  circle  indicates  the  position  of  the  measurement  circle.  The  lower  panels  show  the  measured  phase  of  the 
y  component  of  the  electric  field.  The  measured  phase  is  indicated  by  both  colour  and  vertical  displacement.  The  OAM  mode 
is  found  by  counting  the  number  of  branch  cuts  from  —180°  to  +180°,  i.e.  for  |Z|  =  0,  1,  and  2  we  have  0,  1,  and  2  branch  cuts, 
respectively.  The  orientation  of  the  phase  slope  gives  the  sign  of  the  measured  OAM  mode.  Declining  phase  values  when  the 
xy  plane  is  traversed  in  a  right-handed  sense  indicate  a  negative  mode,  and  increasing  values  a  positive  mode.  The  distortions 
visible  in  the  measured  phase  are  due  to  increased  reflection  when  the  receiving  antenna  is  close  to  the  floor.  The  color  bar  at 
right  shows  phase  in  degrees;  both  the  upper  and  lower  panels  use  the  same  color  scale. 


mode,  as  can  be  expected  from  the  variation  of  the  beam 
maximum  around  the  beam  axis  in  the  radiation  pattern 
for  l  =  1  and  2  in  fig.  2. 

Fig.  4  shows  the  measured  l  spectra,  computed  by 
means  of  a  spatial  Fourier  transform  about  the  z  axis  (see 
Eq.  4  above).  Only  minor  errors  arise  at  the  transmitting 
side,  where  the  feeding  network  delivered  a  maximum 
phase  error  of  3°  and  a  maximum  amplitude  difference 
of  0.03  dB  to  the  six  transmitting  antennas.  These  anten¬ 


nas  were  well-tuned  and  less  than  2%  of  the  power  was 
reflected  back  to  the  transmitter.  The  spread  in  the  spec¬ 
trum  is  primarily  due  to  the  finite  number  of  transmit 
antennas  and  the  reflections  in  the  measurement  cham¬ 
ber.  The  spectra  confirm  that  the  intended  OAM  modes 
were  transmitted  and  correctly  detected. 

In  summary,  we  have  generated  several  radio  OAM 
modes  using  a  circular  antenna  array  and  successfully 
verified  these  modes  via  measurements  of  a  single  elec- 


FIG.  4.  (Colour  on-line)  Angular  momentum  spectra.  Five 
data  sets  of  90  phase  measurements  each  were  taken  every  4° 
around  the  measurement  circle.  Each  set  was  Fourier  trans¬ 
formed,  from  ip  space  to  l  space.  The  magnitudes  of  the 
Fourier  components  corresponding  to  OAM  mode  numbers 
—4  <  l  <  4  are  shown.  The  spectra  show  some  leakage  be¬ 
tween  mode  numbers  but  the  peaks  are  well  correlated  with 
the  intended  OAM  modes,  —2  <  l  <  2.  The  magnitude  scale 
is  normalized  for  each  mode  such  that  the  sum  of  the  magni¬ 
tudes  from  l  =  —  45  to  l  =  +44  is  equal  to  1. 


trie  field  component  using  an  ordinary  electric  half-wave 
dipole;  the  detection  of  OAM  does  not  require  the  use 
of  a  receiving  array  or  a  measurement  of  the  full  electric 
and  magnetic  field  vectors.  In  addition,  we  have  demon¬ 
strated  that  it  is  sufficient  to  measure  the  phase  rotation 
to  determine  the  OAM  mode  of  a  radio  beam.  Although 
the  power  pattern  for  the  |Z|  =  2  beams  was  not  ideal,  the 
spatial  phase  patterns  are  insensitive  to  imperfections  in 
the  beam,  and  the  measured  phase  depends  only  weakly 
on  azimuthal  field  strength  fluctuations  and  on  the  accu¬ 
racy  of  the  circular  measurement  path  around  the  beam 
axis  [6,  19]. 

This  experiment  opens  up  new  possibilities  for  radio 
communications  and  for  radio  and  radar  remote  sensing 
of  rotational  phenomena  in  the  atmosphere  and  space. 
Backscatter  from  OAM-modulated  radar  observations  of 
the  ionosphere  can  be  analysed  using  the  methods  pre¬ 
sented  here  to  diagnose  possible  rotational  properties  of 
plasma  processes.  Artificially-induced  and  natural  radio 
emissions  from  the  ionospheric  plasma  [13,  25,  26]  could 
be  analyzed  for  possible  OAM  effects,  lending  clues  to 
plasma  processes.  The  large-array  radio  telescopes  cur¬ 
rently  being  designed  and  constructed  could  be  used  to 
search  for  signs  of  electromagnetic  OAM  arriving  from  as- 


4 

trophysical  plasmas[23,  27].  Each  OAM  mode  can  act  as 
an  independent  channel  for  transmission  and  reception, 
suggesting  the  possibility  of  increasing  the  information 
transfer  rate  within  existing  measurement  and  communi¬ 
cations  bands.  In  all  applications,  the  reception  of  weak 
radio  signals  located  close  to  an  undesired  strong  source 
might  be  improved  by  orienting  the  the  central  null  in 
an  OAM  receiving  antenna  towards  the  undesired  source. 
For  example,  radio  observations  of  the  solar  corona  could 
be  performed  by  placing  the  central  null  in  an  OAM  beam 
over  the  radio-bright  disk  of  the  sun. 
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